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Abstract 

Many real networks are equipped with short diameters, high clustering, and power-law degree 
distributions. With preferential attachment and network growth, the model by Barabasi and Albert 
simultaneously reproduces these properties, and geographical versions of growing networks have 
also been analyzed. However, nongrowing networks with intrinsic vertex weights often explain 
these features more plausibly, since not all networks are really growing. We propose a geographical 
nongrowing network model with vertex weights. Edges are assumed to form when a pair of vertices 
are spatially close and/or have large summed weights. Our model generalizes a variety of models 
as well as the original nongeographical counterpart, such as the unit disk graph, the Boolean 
model, and the gravity model, which appear in the contexts of percolation, wire communication, 
mechanical and solid physics, sociology, economy, and marketing. In appropriate configurations, 
our model produces small-world networks with power-law degree distributions. We also discuss the 
relation between geography, power laws in networks, and power laws in general quantities serving 
as vertex weights. 

PACS numbers: 89.75.Hc, 89.75.Da, 89.75.Fb 
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I. INTRODUCTION 



Networks of interacting agents such as humans, computers, animal species, proteins, and 
neurons have been investigated vigorously. They are typically complicated, meaning that 
their structures are far from absolutely regular or entirely random. Two principal quantities 
characterizing networks are the average shortest path length L and the clustering coefficient 
C. The number of edges in the shortest path averaged over all vertex pairs defines L. Most 
real networks have small L, namely L oc logn or even less, where n is the number of vertices. 
The local clustering coefficient is the normalized number of connected triangles containing 
a specific vertex. If the vertex degree is fc, or there are k edges adjacent to the vertex, the 
normalization factor is k{k — l)/2. This quantity averaged over all the vertices defines C, 
and real networks usually have large C . A small L and a large C cannot be simultaneously 
realized either by lattices, trees, or the ordinary random graphs Then, Watts and 

Strogatz proposed the small-world networks that fulfill these requirements at the same time 

a. 

Another important observation is that not all but many real networks have power-law 
degree distributions p{k) oc k~"' , typically with scaling exponent 2 < 7 < 3 ■i]. The small- 
world networks are short of the scale-free property. In light of this, Barabasi and Albert 
(BA) proposed a network model that generates scale- free networks with 7 = 30]. Two 
essential features of the BA model are (i) network growth realized by sequentially adding 
vertices and edges, and (ii) preferential attachment, meaning that newly introduced edges 
are more prone to be linked to vertices with larger k. Since the proposal of the BA model, 
its various extensions and related models, such as the fitness model and the hierarchically 
growing models, have been presented. These models are successful in incorporating more 
realistic aspects of networks including tunable 7 and large C that the original BA model 
actually lacks 

To sum up, some BA-type models and hierarchical networks own large C, small L, and 
scale- free p{k). However, real scale- free networks are not necessarily growing. The number 
of vertices may not change greatly over time for networks of friends, companies, interacting 
proteins, and neurons, to name a few. In view of this, a class of nongrowing scale-free 
networks has been studied in which whether edges are created relies on interaction of vertices 
with intrinsic weights Weights represent the fitness of vertices to win edges 
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12l | and are interpreted as, for example, capitals, social skills, activity levels, information 
contents, concentration or mass of physical or chemical substances, and the vertex degree 
itself. The role of such vertex fitness was argued in growing models as well jsl- To our 
surprise, scale-free p{k) emerges even from weight distributions devoid of power laws jlfll 



surprise, sc 



As a remark, if an edge exists when the sum of two vertices' weights exceeds a 
prescribed threshold, the network is equivalent to the so-called threshold graph jl^. This 
case eases analytical treatments. 

Our focus in this paper is the geographical extension of the nongrowing scale-free net- 
works, which has been overlooked so far. Actually, real networks are often embedded in 
topological spaces. Even the Internet, in which the speed of information transmission is 
technically independent of the physical distance, is subject to geographical constraints be- 
cause of wiring costs Il7| . In addition, it is often advantageous to map nonphysical 
quantities or networks into geographical spaces by, for example, the principal component 
analysis. Then, influence of the distance between graduated traits is questioned. 

In fact, the Watts-Strogatz small-world network already addressed this issue since it 
is constructed on lattice substrates Let us denote by g{r) the probability that two 
vertices with distance r are connected. In lattice networks supplied with additional edges, 
where q(r) oc , generated networks have small L when 6 is smaller than a critical value 
Um OtKe.w.e, ,.oba> con.ecHon. a. .00 sca.ce to eU* s.a.l-worM 
The same is true for growing scale-free networks. Although the BA model is irrelevant to 
embedding spaces, which is actually the main cause for small C, it has been extended to 
incorporate underlying geographical spaces and g{r) = r~^. Then, a transition from the 
scale-free to nonscale-free regime as well as one from small L to large L occurs at a certain 

Although g{r) plays a key role in determining the network structure, characterization of 
g{r) in real-world networks still seems controversial. In applications such as the Internet 
routing 21 1 and neural networks j^], g{r) decaying exponentially or in a Gaussian manner 
is commonly used. Exponential decays are also inferred from biological neural networks 
I23I 1. However, many other data are in favor of g(r) oc r~^. For instance, a recent extensive 
analysis of the Internet concludes g{r) oc j6|. Power laws also hold for macroscopic 
connectivity of brain regions identified by correlated activities [g{r) oc r~^] and for 



microscopic neural networks j2^. In social sciences, evaluating g(r) seems more difficult 
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aw and exponential forms 
In the face of the 



because of presumably larger noises. Accordingly, both power- 
of g{r) have been inferred, sometimes even from identical data 
ambiguity of available data, it is worthwhile to examine models to see how various types of 
g{r) affect network properties to help interpret real data. 

In the context of nongrowin g ge ographical networks, there is an algorithm that generates 
p(k) = k~"' with a prescribed 7 flT]. However, investigations of nongrowing geographical net- 
works are largely missing, particularly when interaction of vertices, which is not considered 
in takes place. We examine a geographical threshold network model with various config- 
urations. In Sec. m we review the nongeographical threshold model with vertex weights. In 
Sec. mil we introduce the geographically extended model and analyze some practical cases, 
including the unit disk graph and the gravity model. Section IIVI is devoted to discussing 
our model in the context of network search problems and real data. 



II. NONGEOGRAPHICAL THRESHOLD NETWORK MODEL 

Before taking geography into account, we briefly summarize the ordinary threshold 
network model, which constitutes a subclass of networks with intrinsic vertex weights 

0Q,Q,li3. 

We prepare n vertices denoted hj Vi {1 < i < n), each of which carries a weight variable 
Wi E M. randomly and independently distributed as specified by the density function f{w). 
As mentioned in Sec. HI Wi quantifies the propensity for Vi to gain edges. Let 

/w 
f{w')dw' (1) 
-00 

be the cumulative distribution function. We explain with additive weights since multi- 
plicative weights are transformed into additive weights by taking the logarithm. In the 
nongeographical threshold model, an edge exists between Vi and Vj {i 7^ j) if and only if 
Wi + Wj > 9. When n is sufficiently large, the weight w uniquely determines the vertex 
degree k by 

k = n[l- F{e -w)]. (2) 
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Using Eq. Q, the degree distribution p{k) {0 < k < n) is written as 

p{k) = f{w)'^ 

^ f[0-F-' (1-^)] 

Because the model is simple, L, C, and the correlation between the degrees of adjacent 
vertices can be analytically derived as well Q,^^. The small- world properties characterized 
by a large C and a small L are fulfilled for a wide choice of f{w). More microscopically, 
vertices with small degrees have C{k) near 1 and form the peripheral part of the network. 
It is connected to the cliquish core with larger k and smaller C{k). Strictly speaking, the 
core consists of the vertices with w > 6/2, and the separability of this kind is known in the 

n ~ n 

graph theory A similar separability is also mentioned in other literature 

The degree distribution depends on f{w). An easily solvable example is the exponential 
weight distribution given by 

f{w) = Xe-^"" (0 < w). (4) 

We set 6* > because otherwise the network becomes the complete graph. Although f{w) 
in Eq. is not reminiscent of the power law, substitution of Eq. (jH) into Eq. Q yields 
p{k) oc A;~^ [10]. It follows that C{k) oc and k{k) oc k~^, where k{k) is a measure of 
degree correlation, namely, the average degree of the neighbors of a reference vertex with 
degree k imUM. The same scaling law is also maintained for the logistic distribution 
f{w) = Pe (l + e ^"') , which is just a slight modification of Eq. Another 
major class of fiw) is the Pareto distribution defined by 

/W = -P)"' (5) 

where a,Wo > 0. Equation (0) leads to p{k) oc k~'^ with 7 = (a + l)/a > 1, C{k) oc k^^, and 
k{k) oc k~^ Particularly, C{k) oc k^^ is more consistent with real data compared with 
C{k) oc k~^ derived from Eq. (@]). The asymptotics is the same for the Cauchy distribution 
fiw) = l/[7r{l + w"^)] {w G M), which is devoid of the lower bound of w. The inverse 
problem to determine f{w) from p{k) has also been addressed jl^ . 

A crux of the threshold model is that scale-free p{k) results from a wide class of fiw). 
Analytical and numerical evidence indicates that 7 = 2 is the baseline scaling exponent of 
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the threshold model, which contrasts with 7 = 8 for the BA model [15[. Since the effect of 
a lower bound of w seems marginal, we mainly use the exponential and Pareto distributions 
for the geographically extended model. 



III. GEOGRAPHICAL THRESHOLD NETWORK MODEL 



To generalize the model introduced in Sec. |n] in the geographical case, we assume that 
vertices are uniformly and independently distributed with density p in a ci-dimensional Eu- 
chdean space whose coordinates are denoted by Xi, X2, . . ., Xd- Then a pair of vertices with 
weights w, w', and Euclidean distance r are connected if and only if 



{w + w')h{r) > e, 



(6) 



where h{r) is assumed to decrease in r, although h{r) increasing in r has also been considered 



in other models |4 , 
Boolean model ]2l 







18l |. As a special case, Eq. (jH)) with h{r) oc r ^ is equivalent to the 



Based on Eq. (jHl), two vertices with weights w and w' are adjacent if 

r<h-'( 

\w + w\ 

For a vertex with weight w, the degree k is represented by 



k 



f (w') dw' 



number of vertices in a ball of radius = h 



-1 



e 



w + w' 



(7) 



(8) 



This recovers a general formulation Q| , in which k is calculated from the joint probability 
as a function of w and w' that a pair of vertices are connected. Combination of Eq. (jH)) and 
f{w) provides p{k). If we take an average over w but not over r, we obtain g{r). Although 
g{r) decreases in r if h{r) does, it generally differs from h{r). 



A. Unit disk graph 

If f{w) = S{wo), where 6 is the delta function, two vertices are adjacent when 2woh{r) > 
9. Since h{r) decreases in r, this condition is equivalent to r < ro where 2woh{ro) = 9. 
Accordingly, 

1 (r < rn), 

9{r) ={ - ° ' (9) 

(r > ro). 



and the generated network is the unit disk graph, which is apphed to modehng broadcast 



and sensor networks [16|,|29|. If f{w) has a finite support, the network resembles the unit 
disk graph in the sense that there exists an upper hmit r = ro only below which g{r) > 
0. With this case included, long-range edges are entirely prevented, and the network has 
L oc n^^^, spoiling the small-world property. However, if we allow g{r) = p {r > ro) with 
0<p = 'ri~^^l, we have a type of the Watts-Strogatz small- world networks with small L 
P]. Even so, p{k) is essentially homogeneous. To introduce the scale-free property, we need 
to use more inhomogeneous vertex weights. 

B. Exponential weight distribution with h{r) oc 

Let us consider the exponential weight distribution given in Eq. Q and set 

h{r) = r-^, (10) 

where /3 > 0. This case generalizes the nongeographical model explained in Sec. m which 
corresponds to (3 = 0. For a larger jS, geographical effects are more manifested. As a function 
of the weight, the degree is calculated as follows: 

/•oo /• 

k{w) = / f{w')dw' / p dxi . . . dxd 

Jo J {w+w')/rl^>e 



P 



cie^-rll + l.Au-), (11) 



where 



oo 

r(a,x) = I e-^e-'dt. (12) 

J X 

is the incomplete Gamma function, 

,.^rf^iV (13, 



and r(a) = r(a;,0) is the ordinary Gamma function. To obtain p{k) from k{w), we just 
need to eliminate w from Eq. (fTT|) as we have done in Eq. Q. 

However, an analytical form of p{k) corresponding to Eq. (jlip is unavailable due to 
the incomplete Gamma function. Accordingly, let us deal with some special cases. By 
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integrating Eq. (fT^ by parts, we derive 

r(a,x) = (a-l)!e--'^^ (a G Z), (14) 

a'=0 

where Z is the set of integers. In the hmit /3 — 0, Eq. (fT^ imphes 



a—l , 
X 



lim ^\ ^ 1- = 1. (15) 



Actually, explodes as /5 — ^ because Eq. (fTHjl means lim^^o T ((c^//?) + 1, Aw) = oo, 
reflecting the density notation of the vertex distribution. Putting aside this nonessential 
point, r {{d/P) + 1, Xw) is asymptotically independent of w, and we have 

k{w) oc e^"" (16) 

and 

p{k) oc e-^^"" oc k-^, (17) 
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13, hj. For a suffi- 



which reproduces the results for the nongeographical counterpart 
ciently small (3, Eq. (fT3j) effectively approximates the incomplete Gamma function. Conse- 
quently, scale-free p{k) with 7 = 2 or a slightly larger 7 is almost preserved. 
When (3 = d, we obtain 

k{w) = Cie^'^T {2, Xw) = ci{l + Xw) (18) 

and 



m ^ ^ (19) 

The degree distribution now has an exponential tail, and hubs are less likely compared with 
Eq. ()17|1 . Another special case with (3 = d/2 leads to 

k{w) = ci (2 + 2Xw + A^w^) (20) 

and 

\^-^w exp ^1 — -^^^^j 

"^^^ = 2c,{X + X-w) = 2crVk^i ■ ^''^ 
Equation ()21|) is a stretched exponential distribution with a minor modification factor /c"^/^, 
and p(A;) decays more slowly than in Eq. (fT^ naturally because /3 = d/2 < d. 



Similarly, k and w are connected by a power law relation if /3 > 0. Then, p{k) is a type of 
stretched exponential. In geographical preferential attachment models, the crossover from a 
power-law tail to a stretched exponential tail occurs at a finite value of the control parameter 
similar to P jj, 0, 0, H]- We could say that, in our model, the same transition happens at 
f3 = 0. However, the gist is that for a sufficiently small (3, p{k) is practically indiscernible 
from the scale-free distribution. 

Since it seems difficult to analytically calculate other network characteristics such as L 
and C, we resort to numerics. We uniformly scatter n = 10000 vertices in a two-dimensional 
square lattice with side length 100 and periodic boundary conditions. Because more edges 
obviously means smaller L, the mean degree denoted hjjk) is kept at 20. The analytic 
expression for (k) is available only when /5 = as follows |l3L Il5l|: 

{k) = e-^' {pi'' + Xe) , (22) 

where / is the side length of the area. Therefore, we manually modulate 6 to preserve (k) 
as we vary j3. Excluding isolated components, which actually consist of just a few vertices, 
we show a dependence of L and C on /5 in Figs. Ufa) andHJ^b), respectively. Although the 
main simulations are done with n = 10000 (thickest lines), we also show results for n = 2000 
(thinnest lines), 4000, 6000, and 8000. The inset of Fig. ^^a) shows the dependence of 
L on n, with upper lines corresponding to larger values of (3. Figure ^a) shows that L 
is insensitive to n only when j3 < 0.5. We expect that L oc \ogn approximately holds in 
this regime. On the other hand, we expect L oc n^^'^ or similar scaling for larger p. As 
f3 increases, C decrease to some extent but not too much to spoil the clustering property 
[Fig. Cl^b)]. Figures d^c), [H^d), andlH^e) show p(A;) (crosses) and C{k) (circles) for (3 = 0.5, 
P = 1.5, and f3 = 2.5, respectively. As expected, small (3 yields a long tail indicative of the 
power law [Fig. C^c)]. In contrast. Fig. ^e) shows that p{k) decays much faster when jS is 
larger. Consequently, networks generated by sufficiently small (3 are endowed with the scale- 
free and small- world properties simultaneously in a geographical context, which extrapolates 
the nongeographical results with = 0. In regard to the vertex- wise clustering coefficients, 
C{k) oc k~'^ holds when /? = The numerical results in Figs. [U^c), [H^d), andd^e) 

(circles) support C{k) oc except that vertices with small C{k) are more scarce for larger 
p. 
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The probability g{r) that two vertices with distance r are adjacent becomes 
g{r) = 



(w+w')/rP>e 



Xe-^'^dw I Xe-'''"'dw' + I Xe-^'"dw 
e->^s-\xer^ + 1), 



(23) 



indicating a stretched exponential decay in r unless /? = 0. Particularly, the main decay 
rates for (3=1 and /? = 2, respectively, correspond to the standard exponential and the 
Gaussian, which are widely used in applications 21, 2^. As a general remark, g{r) does not 
coincide with h{r) oc even asymptotically. 

The loss of the small-world property for large (3 seems to stem from the (stretched) expo- 
nential deca y o f g{r). In addition, g{r) derived here qualitatively disagrees with many real 



data 
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l2q . As a result, exponential types of g{r) and the Gaussian g{r) may be far 
from universal. This is a striking caveat to many fields, such as neuroscience, social dynam- 
ics, and epidemics, which conventionally assume geographical networks with exponentially 
decaying or Gaussian g{r). We do not explore consequences of h{r) that decays faster than 
h{r) oc r~^, since such an h{r) must yield even larger L. On the other hand, h{r) with 
slower decays, or h{r) oc (logr)~^, is examined in Sec. IIII Dl 



C. Power-law weight distribution with h{r) oc r 



Quantities that can serve as vertex weights, such as the ci ty a nd firm sizes |3Cil3lLl32ll33|. 



number of pages in a website j3^, land prices j35|, incomes 36 1, importance of airports 9|, 
and importance of academic authors , are often distributed according to power laws. The 
history of these power laws is much longer, dating back to the Pareto and Zipf laws, than 
those recently found for networks j^. The simplest way to associate the power laws of 
networks with those of vertex weights is simply to interpret the vertex degree as the weight. 
However, w and k are generally nonidentical P, Q| . 

Let f{w) be the Pareto distribution given in Eq. (0). With the interaction strength 
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decaying algebraically [Eq. (fTUj)]. we have 

' wo/w 

Convergence of k{w) necessitates —a + d/P < 0. In the limit w — oo, it holds that 



-r f ^ + l) ^('^Z/^)- r ^^^S' dx {w>wo). (24) 



a 

W ' 



Therefore, 



and 



Il±^rfa; oc r -^dx oc ( — ) . (25) 

Wo/w J Wo/w / 

k{w) oc w'^/^ (26) 

_a_ (wo)"-'^^ 



In contrast to the stretched exponential scenario clarified in Sec. IIII B| the power-law weight 
distribution produces scale- free p{k) = with 7 = 1 + (aP/d). 
For r large enough to satisfy 6r^ > 2wq, 



+ / x-'^-^(6-x)-'^(ia;, (28) 



— Wo, 

where h = 9r^/wQ. To show that the integral in Eq. pH|l tends to be proportional to r~"^ 
as r — i> 00, let us evaluate b"" ^ A{x)dx, where A{x) = x^"'^^{b — x)^". First, we obtain 

liminffe" / A( > lim / x'^'-^dx 

= - lim h-ib-iy] =-. (29) 
To bound the integral from the above in the limit b —* 00, let us assume b > 4. Noting that 
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A{x) takes the minimum at x = (a + l)6/(2a + 1) and that (PA{x)/dx'^ > 0, we derive 
A{x)dx 



Hm sup 6" 

fe— too 



< hm 





' 36 " 


y(fe+2)/3 


[2(6-1) J 



X 



6" 

-"-^dx + — 



A 



+- 



hm 



A 



(a + 1)6 
2a + 1 
36 



+ 



2(6-1 
6 



+ A(6-1) 
3 



(6-1) 



6 + 2 
3 

(a + 1)6 
2a + 1 



;a + l)6 
2a + 1 



1 



6 + 2 



2 \b-l 
1 /3^ ° 
a I 2 



a6 — 2a — 1 



+ 



_(2a+ 1)(6- 1) 
a 

< 00. 



O (6-") 



(a + 1)6 6 + 2" 



2a + 1 



2(2a + 1 

Equations and (jHUj) ahow us to conclude an algebraic decay g{r) oc r~"^ in contrast to 



(30) 



Eq. Discussion on network structure is postponed to Sec. IlllPj where we analyze the 

gravity model, which ends up with the same asymptotic behavior of p{k) and g{r). 



D. Gravity model with Pareto f{w) 

As shown in Sec. IIIIB| given the exponentially distributed w, h{r) oc with a suffi- 
ciently small /3 yields more or less desired network properties. More rapidly decaying h{r) 
makes g{r) decrease too fast to elicit small L. How about h{r) that decays more slowly? 
To address this issue, we apply h{r) oc (logr)~^ Since logr can be negative, let us rewrite 
Eq. (0) as 

w + w'>9\ogr. (31) 

Equation (p?T| is equivalent to 

e'^e'"' > (32) 

Since edge formation is suppressed by increasing either P or 9, let us reinterpret 9 in Eq. 
as j3, which does not essentially change the model. Further rescaling of the parameters by 
W = e"^, W = e"^' , and R = 9-^/l^r transforms Eq. ^ into 

WW 



> 9. (33) 



This is the gravity model often used in physics, sociology, economics, and marketing 
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38, 



26 



30 



The gravity model is suitable in describing interaction of particles in geographical 

12 



spaces when the physical gravity (/? = 2) or similar mass interaction based on, for example, 
populations or chemical substances, is active. In the sociological context, the original model 
stipulates 3 = 1 |8fl| . but (3 ranging from 0.2 to 2.7 have been inferred from later real data 



:es 3 = 1 

^, H H Isa i4i 



The original gravity model is geographical but neglects weight distributions. On the 
other hand, multiplicatively interacting weights with power-law f{w) are used to generate 
solvable scale- free networks, but they ignore geography 0| . We are interested in combined 
effects of geography and dispersed vertex weights. The transformation from Eq. (jHTjl to 
Eq. (jnni) also rescales f{w) unless it is the delta function. When the weights in Eq. (p?T|l 
follows the exponential distribution given in Eq. (^, the density f{W) of the weights in 
Eq. (jnSl) becomes 

dw , / 1 ^ 



nw) = fH- = X^-) , (34) 

namely the Pareto distribution with a = X and wq = 1- Although we have started with 
h{r) oc (logr)~^ and additive weights, we switch to the gravity-model notation for conve- 
nience. Now we rewrite Eq. fl33|) as 

^ S " (3«) 

and investigate the network structure when f{w) is the Pareto distribution. 

Before moving on to the Pareto f{w), let us note that f{w) with a finite support only 
allows local interaction, as explained in Sec. IIII Al Then the gravity model yields L oc n^^^, 
which is realized by atomic and molecular interaction by centrifugal or electric forces; they 
practically interact only with others nearby. With the Pareto f{w), which facilitates more 
global interaction, we obtain 



dw' 



I wo Wo \W' 

C2 w'^lf', (36) 



C2 = ^577^- — j< r - + 1 . (37) 



where 

Equation (jHBj) is essentially the same as Eq. (j211), and [3 > d/a must be satisfied for C2 > 0. 
The original gravity model for social interaction has {3 = 1 and d = 2 [s^ , and hence a > 2 
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is necessary. Combination of Eqs. (0) and (IHUj) yields 




Oj^-l-ia/3/d) 



a 



(38) 



When > w^/O, we obtain 
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(39) 



Comparison of Eqs. ^T7\\ and (PHj) with Eqs. (jHHj) and (jH^ reveals that the asymptotic 
behavior of p{k) and that of g{r) coincide with those of the additive weight model with the 
Pareto f{w) and h{r) = r^^. Given the Pareto f{w) and h{r) = r~^, whether multiplicative 
or additive interaction is used does not matter so much. 

Numerically evaluated L, C, p{k), and C{k) for varying /? are shown in Fig. El We set 
n = 10000, a = 1, t^o = 1, and 



constant at 20. Figures Efa) andlSfb) show that L and C have a similar dependence on (3 to 
the additive weight model with exponential f{w) [Figs.IH^a) andHJ^b)]. FigureEfa) indicates 
that a transition from a small-L regime to a large-L regime occurs somewhere around (3 = 3. 
Since Fig. I2tb) supports that C remains finite for large n irrespective of (3, the small-world 
property is suggested for small (3. The transition appears similar to the phase transition in 
geographical BA models 0,0]. However, in those models, 7 does not change in /3 > as far 
as the network is in the small- world regime, whereas it does change here (but see Q]). As 
shown in Figs. Efc), El^d), and|21^e), p{k) (crosses) obey power laws whose scaling exponents 
are well predicted by Eq. (jHH)) (lines). Consequently, the weighted gravity model reahzes 
scale-free small-world networks when f3 is small enough. In this scheme, 7 is tunable by 
varying a, (3, and d. Circles in Figs. |2tc),|2fd), and|2fe) indicate C{k) oc k~'^' with 7' = 2 or 
somewhat smaller. Finally, numerically obtained g{r) shown by circles in Figs. Efa) {(3 = 2) 
andOfb) (/3 = 3) decays algebraically as predicted by Eq. (jHEI)- 

A generated network is shown in Fig. |3]for n = 100, d = 1, a = 1, wq = 1, (3 = 1, and 
hence 7 = 2. For demonstration purposes, the vertices are aligned on a one-dimensional 
ring. In spite of the small size, the figure is indicative of the scale-free and small-world 




(40) 
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1-world networks on a 



3 



properties. It is visually comparable to the BA-type scale-free sma^ 

n 

ring [5j and the Watts-Strogatz non- scale-free small-world networks 

In other geographical network models, L becomes large if g{r) decays faster than g{r) oc 
with a certain 5 > 0. For example, a nonscale-free weightless network model on a 
lattice owns an ultrasmall L = 0(1) for 6 < d, small L = O(logn) for d < 6 < 2d, and 
large L = 0{n^^'^) for 6 > 2d In another nonscale-free network, 6 = d + 1 divides the 
small- world and large- world regimes Also in a one-dimensional geographical scale- free 
network model with preferential attachment, a similar phase transition occurs at 5 = 1 P,!?!- 
Based on Figs, ^a) andlSfa), we anticipate that the gravity model has the phase transition 
at a critical value 6 under which the network is geographical, scale-free, and small-world 
at the same time. We do not examine h{r) decaying faster than (logr)~^ in the additive 
weight notation [Eq. (jHH)] or equivalently h{r) decaying faster than algebraically in the 
multiplicative weight notation [Eq. ()35|)]. for which we expect too large L. Let us mention 
that h{r) oc (logr)~^, which other models have largely neglected, may be appropriate if 
weight interaction is effectively additive. 

The results in Sec. IIII (^l and those in this section can be captured as a spatial extension 



of the results in |lj], which addresses the inverse problem to determine /{w) from p{k). To 
obtain p{k) oc fc"'^, a pair of vertices with weights w and w' that follow f{w) = Ae~^"' with 
A = 1 are connected with probability proportional to exp [—{w + 7 + 1)] 14j. In the 

gravity model, we have defined W = exp{w) and W = exp(w') so that W and W follow 
the Pareto distribution. The probability that the two vertices are connected is proportional 
to the volume of a d- dimensional ball with radius ro, where WW/tq = 9. This probability 
is proportional to Tq oc (WW /O)'^^'^ oc exp[{w + w')d/(3]. We should have d/l3 = 1/(7 — 1), 
which is consistent with Eq. fl38|) since a = A = 1. 



E. Gravity model with exponential f{w) 

Let us treat the gravity model with the exponential weight distribution. This configu- 
ration is equivalent to the model with additive weight interaction, h{r) oc (logr)~^, and a 
weight distribution less broad than the exponential distribution. It follows that 

k{w) =pj^ Ae-^^'V'^/^r + 1^ dw' = CiF + 1^ w"/^ (41) 
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and 



p{k) 



(3X 



exp 



-A 



k 



ciT {d/(3 + 



13/d 



(42) 



which is a stretched exponential with a modifying factor k i+^Z*^. Furthermore, we have 



cos 



(m2 + 4A2)3/2 



du 



AX' 



dt = AX'^Ko{2XVe^), 



(43) 



where -^'0(2;) is the modified Bessel function of the second kind j4l', pp.185, 187-188, and 
206]. Since Kq{x) tends to 



TT 



Ko(x) = W — e 




l-:- + - 



6X 



as X 



00 



41 



p. 202], Eq. pnj) asymptotically behaves as 



00). 



(44) 



(45) 



With the arguments in Sec. HIT Dl taken into account, Eq. ()45p implies that g{r) decays too 
fast to make the network small- world. A lesson is that f{w) considerably influences network 
properties, which is not the case for the nongeographical counterpart 15[- Particularly, the 
Pareto f{w) can yield scale-free p{k) and the small-world properties, regardless of whether 
weight interaction is additive or multiplicative. On the other hand, the exponential f{w) 
explored in this section and Sec. IIII Bl induces exponential types of p{k) and large L. 



IV. SCALE-FREE NETWORKS AND SCALE-FREE WEIGHT DISTRIBUTIONS 

Among the configurations considered in Sec. lIIH the additive weight model and the gravity 
model with scale- free f{w) and scale- free h{r) generate small- world networks with scale- free 
p{k). In this regime, our model relates scale-free p{k), which is of recent research interest, 
to general power law distributions in nature that have a long history tracing back to Pareto. 
Let us discuss the relevance of our model to real data. 

There is a body of evidence that quantities potentially serving as vertex weights are 
distributed according to power laws f{w) oc w~°'~^. For example, the celebrated Pareto and 



Zipf laws dictate that incomes and city sizes follow power laws with a + 1 = 2.0 



30, m^. 



More recent data analyses confirm power laws in countrywise city sizes {a + 1 = 1.81-2.96 
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with mean a + 1 = 2.136) |32j |. firm sizes (a + 1 = 2) |33| . the number of pages per website 
(a + 1 = 1.65-1.91) 34], land prices (a + 1 = 2.1-2.76) 35|, incomes (a + 1 = 1.7-2.4) 
and importance of airports (a + 1 = 1.67) j^, to name a few. On the other hand, the original 
gravity model disregarding weight distributions assumes /3 = 1 and d = 2 j^^. Application 
of the values of a mentioned above to the weighted gravity model yields 'j = l + a(3/d = 1.32- 
1.98, which is too small to fit real network data whose 7 mostly falls between 2 and 3 P|. 
As another indication, an extensive data analysis of the Internet revealed g{r) oc with 

6 = 1 p]. If our model could underlie the Internet, it should mean a(3 = 6 = 1, and hence 

7 = 1 + af3/d = 3/2 since d = 2. This 7 is again too small for the real Internet and related 
computer networks that have 7 = 1.9-2.8 

However, we regard that our model is not necessarily implausible. First, our model and 
also the nonspatial threshold model do not aim to describe growing networks; the Internet is 
a typical example of growing networks [3] • Our goal is rather to discuss nongrowing networks 
in a geographical context. As a supporting example, connectivity networks of brain regions 
have 7 = 2, 5 = a/3 = 2, and d = 2 [2J|, which are roughly consistent with Eq. (jHEI)- 
Actually, the brain network does not grow so much once an animal is born. 

Second, estimation of (3 involves much fluctuation because of the difficulty in data acqui- 
sition. Since the pro pos al of the gravity model in which (3 = 1 was inferred from railway and 
highway travel data |3Q;], analyses of various social activity data have offered a wide range of 
(3. Among them are investigations of air travelsL/3 = 0.2-2.0) HHQ, 

journey to work 

0.5-1.2) 0, migration {{3 = 1.59,2.49) H, cedar rapids direct contacts {{3 = 2.74) 
marriage (/? = 0.59,1.53,1.59) [s^, and memorizable social interaction {(3 = 2) j^]. 
The ambiguity and the activity dependence of (3 render the evaluation of 7 pretty uncertain. 
Precision of (3 in classical studies was also low because of small data sizes. To undertake 
more detailed and large-scale data analysis as in j6, 40] is important. 

Third, the interaction strength, which is assumed to be proportional to wvwo/r^ in the 
gravity model, may be nonlinear in weights. For example, use of lyfiyf/r^ j2a] results in 
7 = 1 -|- aP/xd. Real data actually support 0.73 < x < 1.05 ^^l, and x smaller than 
1 increases 7 to make it more realistic. By the same token, replacing {wi + W2)/r^ with 
{wi + W2Y Ir^ in the additive notation effectively changes (3 to 13/x. 

Next, let us relate our model to network search problems in which an agent on a vertex 
attempts to reach an unknown destination by traveling on edges. In small-world networks 
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defined by lattices supplied with long-range connections with density g{r) oc r^^, which are 
essentially equivalent to the random connection model [2^, optimal search performance is 
realized when 6 = d Even though the weighted gravity model is a different model, 

simple adoption of our formula suggests d = 5 = a(3 and = 1 + a[3 / d = 2. Computer- 
related networks usually have 7 > 2 presumably because they are growing. However, some 
social networks and peer-to-peer networks, which may be considered to be nongrowing, own 
7 close to 2 enhancing the search ability. 

Similarly, emergence of small- world networks in a geographical framework requires d+1 > 
6 > d, while latticelike networks result from 6 > d+ 1 , and 6 < d induces randomlike networks 
with low clustering jig] . Simple-minded substitution of 5 = a/3 leads to d + 1 > ajS > d 
and 2 < 7 < 2 + d~^. Since we usually have c? = 2 or 3, 7 associated with general 
nongrowing small-world networks may be close to 2. To summarize, scale-free networks 
with 7 around 2 may be optimal in the sense of the search performance and the small- 
world property. In addition, 7 = 2 is the baseline scaling exponent of the nongeographical 
threshold g raph J1 5 t] , and it may also be the case for general cooperative nongrowing networks 
Ifll U, 12, llJ. Ilfil Q|. In contrast, 7 = 3 is an important phase-transition point for 
percolation and dynamic epidemic processes |42|. The BA model has 7 = 3, which may 
set the baseline 7 for other competitive growing networks . Our current 

speculation stems from the ansatz 7 = 1 + a/S/d = 1 + 6/d plugged into the results obtained 
from other models. Further investigation of this issue is an important future problem. 



V. CONCLUSIONS 



We have proposed and analyzed a geographical nongrowing network model based on 
thresholding the sum of two vertex weights. Our model contrasts with geographical growing 
models based on the BA model, and it naturally extends the threshold graph, the unit disk 
graph, and the gravity model, which are widely used in a range of fields. In proper regimes, 
small-world networks with scale-free degree distributions and the connection probability 
algebraically decaying in distance are generated, and they are consistent with many real 
data. In contrast to the nongeographical threshold model, what weight distribution is used 
matters for network properties. For scale-free networks to emerge, the weight should be 
distributed as specified by power laws. The weight distribution and the degree distribution 
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generally have different scaling exponents, and they are bridged by a relation involving the 
spatial dimension and the decay rate of the interaction strength. 
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Figure captions 

Figure 1: Network properties with h{r) = and the exponential weight distribution 
with A = 1 and (k) = 20. Dependence of (a) L and (b) C on (3 for n = 2000 (thinnest 
hues), 4000, 6000, 8000, and 10000 (thickest hues) is presented. The relation between L and 
n is shown in the inset of (a), with upper lines corresponding to larger (3. Also shown are 
numerically obtained p{k) (crosses) and C{k) (circles) with n = 10000 for (c) (3 = 0.5, (d) 
p = 1.5, and (e) p = 2.5. 

Figure 2: Network properties for the gravity model with the Pareto weight distribution 
with a = 1, Wq = 1, and {k) = 20. Dependence of (a) L and (b) C on /3 is presented. 
Also shown are numerically obtained p{k) (crosses), C{k) (circles) with n = 10000, and the 
theoretical prediction p(k) oc k~^~"-f^/^ (lines) for (c) /3 = 2, (d) /? = 3, and (e) /3 = 4. 

Figure 3: Numerically obtained g{r) (circles) and the prediction by Eq. ()39|1 (lines) for 
the gravity model with (a) /5 = 2 and (b) (3 = 3. The other parameter values are the same 
as those used in Fig. El 

Figure 4: An example of the weighted gravity model on a one- dimensional ring. We set 
n = 100, (3 = 1, and {k) = 6. The Pareto weight distribution with a = 1 and t^o = 1 is used. 
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FIG. 4 
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